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WAVELENGTH  SCALING  FOR  REACTOR-SIZE  LASER  FUSION  TARGETS 


In  order  to  achieve  high  pellet  gains  for  laser  fusion,  one  must  control 

a  number  of  distinct  critical  physics  elements.^  Two  of  these  elements  are 

a  sufficient  coupling  efficiency  of  laser  energy  into  the  Imploding  DT  fuel 

(5-lOZ),  and  a  low  fuel  preheat  (2-4  times  Fermi  temperature).  To  satisfy 

these  conditions  one  must  limit  the  deleterious  plasma  instabilities,  such. 

as  stimulated  backscatter  and  suprathermal  electron  generation.  It  has 

been  widely  suggested  that  this  control  be  accomplished  first  by  use  of  thin— 
2  3 

shell  or  double-shell  pellet  designs  that  use  a  lower  peak  laser  irradlance, 

and  second  by  the  use  of  shorter  laser  wavelengths^  in  the  range  to  %  micron 

2 

since  the  strength  of  many  of  the  plasma  instabilities  scales  as  I  A  . 

This  article  will  present  the  scaling  law  of  another  critical  element: 

the  symmetry  of  the  pellet  Implosion.  For  high  gain,  the  asymmetry  of 

the  ablation  pressure  must  be  less  than>a  few  percent;^  i.e.  6P/P  ^  1-3Z. 

Either  the  laser  illumination  on  the  pellet  must  also  be  uniform  to  this 

percentage,  or  one  needs  some  physics  mechanism  to  smooth  out  the  laser- 

nonunlformlcies.  Nuckolls  et  al.^  proposed  such  a  physics  mechanism:  by 

designing  their  shaped-pulse  pellet  so  that  the  distance  D  between  the 

laser  absorption  surface  and  the  ablation  surface  was  Itself  on  the  order 

of  the  pellet  radius,  (D<vR),  lateral  thermal  conduction  could  smooth  out 

the  laser  nonuniformities,  yielding  dP/P  «dX/l.  Such  smoothing  has  been 

observed  experimentally.^  Using  a  hydrodynamic  computer  code,  we  have  found 

a  scaling  law  for  this  separation  distance  D.  This  D  scales  much  more 

strongly  with  wavelength  than  do  the  plasma  instability  effects.  We  find 
3  8  0  7 

Da/(i  a  ’  )  *  for  a  given  pellet  size.  This  places  a  strong  premium  on 
the  use  of  longer  laser  wavelengths,  in  the  near-infrared,  for  reactor-sized 
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(1-5  mm)  pellets*  Plasma  Instabilities  may  then  have  to  be  controlled  by 
ocher  techniques,  such  as  the  use  of  a  very  broad-band  laser,  and/or  the 
use  of  thinner  pellet  shells.  (Symmetrizatlon  by  the  use  of  x-ray  coupling 
Is  not  considered  in  this  article.) 

Because  the  code  compares  well  with  experimental  data  at  the  Naval 
Research  Laboratory  on  ablation  pressures,  mass  ablation  rates,  and  target 
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velocities  from  the  ablative  acceleration  of  flat  targets,  ve  have 
a  degree  of  confidence  in  our  scaling  result.  Our  code  also  has  some 
distinctive  capabilities,  discussed  below.  We  have  also  applied  the  code 
to  a  study  of  the  higher  density  region  of  the  pellet  when  driven  by 
moderate  laser  intensities,  (^10  W/cm  )  where  fewer  uncertain  plasma 
effects  exist. 

The  Model.  We  used  1-D  spherical  hydrodynamics,  with  time-dependent 
equations  for  continuity  of  mass,  momentum,  and  energy,  including  light 
pressure.  The  laser  light  not  absorbed  by  inverse  bremestrahlung  was 
deposited  smoothly  over  two  zones  including  the  critical  surface,  so  as 
not  to  generate  artificial  shocks. 

Thermal  conduction  was  computed  using  the  usual  smoothing  technique 
with  the  conduction  coefficient  given  by  1/K  »  ^^^l***  from 

Spitzer-Harm,  and  «  f  /kT/nP  kT/|  VI  {  •  The  flux-limiter  is  a 
construct  used  in  fluid  codes,  for  which  there  are  not  tests  yet  on  reactor^ 
sized  plasmas.  We  tried  values  for  f  ranging  front  0.6  to  0.03.  We  have 

I 

presented  results  using  the  Bell-Evans-Nlcholas  value  of  f»0.1  which  ve 

1 

believe  to  be  on  the  surest  theoretical  grounds.  The  separation  distance 
D  is  rather  insensitive  to  f,  for  £>  0.1,  as  long  as  the  Mach  number  at  the 
crlticel  surface  is  also  not  too  much  greater  than  one  (in  other  words, 
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small  spherical  divergence).  For  £  »  0.03  the  flow  Mach  iiuiid>er  Is  llsltcd 
Co  less  than  0.5  in  Che  overdense  region,  and  then  the  separation  distance 
D  is  reduced  by  a  factor  of  3  to  4. 

The  equations  were  integrated  in  strict  conservation  form,  guaranteeing 
correct  motion  of  shock  waves  in  the  accelerated  dense  shell.  Pressure 
and  temperature  were  thus  derived  by  subtracting  directed  kinetic  energy 
from  total  energy.  The  fluid  equations  were  closed  with  equations  of  state 
that  express  P  -  F(E,p)  and  T  »  G(£,p).  Near  the  ablation  surface  these 
equations  of  state  must  include  the  effects  of  solid  material,  Fermi^ 
temperature,  and  fractional  ionization.  Ue  based  our  equation  of  state 
on  the  solid  material  models  used  in  the  analytic  model  (ANEOS)  in  the 
CHARTD  code,^^  and  expansions  of  the  Fermi-Dirac  functions  by  Clayton, 
incorporating  a  phenomenological  cold  compression  curve  to  model  the 
bulk  modulus  and  the  binding  energy  of  the  atoms.  This  model  goes  smoothly 
from  a  Thoznas-Fermi-Dlrac  model  for  degenerate  electrons  at  high  density 
and  low  temperatures  to  a  Saha  equilibrium  at  high  temperatures  and  low 
densities.  Such  a  treatment  is  necessary  in  order  to  correctly  calculate 
the  acceleration  of  near-solid,  low-temperature  target  shells  without  having 
the  shell  compress  or  explode  unrealistically. 

The  equations  were  integrated  numerically  with  a  sliding  zone,  Eulerlan, 
flux-corrected-transport  (FCT)  code  and  solved  with  a  flux<*conservlng 
scheme  such  that  the  fluxes  between  zones  result  from  the  difference  in  fluid 
and  grid  velocities.  The  outer  boundary  zones  moved  in  a  Lagrangian 
manner,  while  the  intermediary  zones  moved  in  a  prescribed  smooth  continuous 
manner.  Near  the  ablation  surface,  the  mass  per  unit  zone  was  chosen 
proportional  to  the  density  (constant  zone  size) •  Beyond  the  ablation 
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surface^  the  mass  per  zone  continually  decreased  until  it  reached  a 
minimum  value.  Thereafter  the  mass  per  zone  was  held  fixed  up  to  the 
critical  surface  (increasing  zone  size).  Outside  the  critical  surface, 
the  mass  per  zone  increased  exponentially.  This  technique  allows  a. 
grid  resolution  of  a  few  tenths  of  a  micron  near  the  ablation  surface; 
it  also  has  the  distinct  advantage  over  purely  Lagrangian  codes  of  not 
wasting  zones  where  not  needed. 

Rather  than  attempt  to  design  a  set  of  optimized  implosions^  we 

chose  pure  CU  shells  of  constant  thickness  of  125  ym.  We  used  no  prepulses 

to  generate  an  artificially  large  initial  separation  distance  D„  since 

so  far  such  prepulses  have  led  to  shock  heating  that  Is  unacceptable  for 

high  gain  pellets.  The  results  here  are  expressed  in  terms  of  absorbed 

Irradlance,  since  the  absorption  process  is  not  well  understood  for 

rc^ctor-slzed  pellets.  (Based  on  current  knowledge,  ve  guess  that  net 

absorptions  will  lie  in  the  range  of  40%  90%  for  the  intensity-wavelength 

regimes  treated  in  this  article.)  Reactor-sized  pellets  will  have  radii 

of  1-5  mm.  Since  laser  energy  scales  as  R  (aR/R)  ,  a  much  larger  radius 

would  Imply  an  excessive  aspect  ratio  or  excessive  laser  energy*  A 

smaller  radius  cannot  contain  sufficient  fuel  for  efficient  burn  and  high 

gain.  The  intensity  profile  was  ramped  linearly  for  2  nsec,  followed  by 
2 

I^(t)R  (t)  •  constant.  The  last  variation  had  no  appreciable  effect, 
since  we  found  our  qtiasi-steady  profiles  by  looking  early  in  time  before 
the  bulk  of  the  shell  had  moved  Inward  appreciably. 

The  Results.  Calculations  at  laser  wavelengths  of  0.53  pm,  1.06 
and  2.7  ym  with  pellet  radii  of  Imra,  2mm,  and  5aun  were  used  to  derive 
scaling  laws  for  the  pressure  P  and  the  separation  ratio  D/R  versus  laser 
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intensity,  laser  wavelength,  and  pellet  radius.  A  least  squares  fit 


of  the  form  P«»  (I^  ^-0.25jj-0. 18^0.7  shown  in  Figure  1.  Figure  2  shows  _ Q 


D/R  versus  laser  Intensity  for  three  wavelengths  and  two  initial  radii. 


He  found  a  least  squares  fit  of  the  form 
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In  the  region  of  interest ,  B  Cf  Immg  and  is  only  weakly  dependent  on 
This  form  goes  to  the  correct  limits  for  both  large  and  small  radii* 
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This  surprisingly  strong  wavelength  dependence  is  the  result  of  two  factors: 

longer  wavelength  light  is  absorbed  at  lover  densities »  and  these  lower 

densities  have  longer  scale  lengths. 

12 

Emery  et  al.  have  shown,  for  flat  geometry,  that  if  the  separation 

distance  D  is  more  than  0.3  times  the  Ibser  inhomogeneity  wavelength,  then 

dP/PaO.l  5I/I*  We  expect  roughly  the  same  factor  of  ten  smoothing  in 

spherical  geometry,  if  D/R  »  0*3.  From  Figure  2,  this  implies  a  minimum 

13  2 

laser  intensity  of  about  0.5  -  2  x  10  W/cm  with  2.7  ym  light,  about 

14  2  IS  2 

0.8  -  2  X  10  W/cm  with  lym  light,  and  more  than  10  W/cm  with  0.5  ym 

light.  Still  shorter  laser  wavelengths,  of  1/3  to  \  micron,  would  require 

much  higher  laser  intensities. 

For  today's  moderate  aspect  ratio  pellet,  with  **  8-20,  an  ablation 

pressure  of  about  10-30  megabars  is  required  to  drive  the  shell  Inward 

to  fusion-level  velocities,  using  a  relatively  unshaped  laser  pulse* 

Extrapolating  from  Figure  1,  we  see  that  the  optimal  laser  intensity  Is 

14  2 

probably  a  few  times  10^  W/cm  •  It  appears  that  a  one  micron  laser  should 
suffice  for  producing  both  sufficient  ablation  pressure  and  sufficient  Z-D 


^  <  tumid 


Fig.  1  ~  Ablation  surface  pressure  scaling,  with  several 
laser  wavelengths  and  radii 
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smoothing.  (For  a  flux  limit  of  £  ■>  0.03  the  optimal  wavelength  would 
increase  to  two  to  three  microns.) 

Discussion 

Why  not  use  a  shorter  laser  wavelength?  To  keep  the  same  smoothing 

factor  in  the  same  size  pellet  with  a  h  micron  laser ^  we  would  need  an  In- 

3  8 

tensity  that  was  increased  by  about  the  factor  2  *  14.  Then  the  plasma 

2 

Instability  parameter  IX  would  increase  by  the  factor  3.5.  Therefore 
shorter  wavelengths,  for  the  same  amount  of  smoothing,  would  actually  increase 
the  potential  for  plasma  instabilities.  Equally  important,  the  ablation 
pressure  at  the  higher  laser  intensity  might  be  too  large,  thereby  driving 
strong  shocks  through  the  shell  that  would  preheat  the  pellet  fuel*  If  one 
reverts  to  a  shaped-pulse  pellet  concept,  one  can  use  these  higher  laser  in¬ 
tensities,  but  only  after  the  pellet  shAl  has  moved  significantly  inward, 
again  with  unacceptable  as^-trste tries. 

What  about  a  much  longer  wavelength  laser?  There  would  be  much  more 

of  a  safety  margin  for  lateral  smoothing  with  a  2-3  micron  laser,  at  the 

14  2 

required  intensity  of  a  few  times  10  W/cm  ,  and  in  fact  this  may  be 

2 

necessary.  Buu  for  fixed  ablation  pressure,  the  instability  factor  I  X 
would  again  increase. 

There  seems  to  be  a  barrier  at  shorter  wavelengths  because  the  smoothing 
requirement  quickly  drives  one  to  unacceptably  high  laser  intensities  and 
ablation  pressures.  And  at  longer  wavelengths,  the  ablation  pressure  require¬ 
ment  on  the  laser  intensity  soon  exceeds  the  threshold  for  plasma  instabilities* 
The  actual  location  of  the  wavelength  window  is  somewhat  uncertain,  but  It 
appears  to  lie  between  ^  micron  and  2.7  microns. 


Since  we  first  proposed  this  wavelength  window  concept »  there  have 

been  several  suggestions  for  widening  the  window.  We  have  suggested  usin^ 

a  very  broad  bandwidth  laser,^  auch  as  HF  at  2,7-3.1  micron,, to  raise 

13 

the  plasma  instability  threshold  by  an  order  of  magnitude.  Fabre  has 

suggested  using  a  much  smaller  pellet  radius  R,  but  it  would  be  difficult 

to  contain  enough  DT  fuel  as  stated  above.  Nuckolls  has  suggested  using 

14 

a  variant  of  his  old  proposal  for  a  series  of  temporal  changes  in  the 
laser  frequency,  starting  at  longer  wavelengths  for  maximum  smoothing,  and 
ending  with  shorter  wavelengths  for  maximum  coupling  efficiency. 

It  is  best  to  finish  with  some  caveats.  First  our  hydrodynamic 
model,  and  any  hydrodynamic  model,  tends  to  fail  when  trying  to  accurately 
treat  plasma  phenomena  such  as  fast  electron  transport  and  filamentation. 

And  the  most  interesting  regime  for  reaotor-slzed  targets  is  Just  at  that 

irradiance  where  hydrodynamic  models  tend  to  fail.  Second,  if  the  long- 

wavelength  laser -plasma  coupling  is  poor,  then  there  may  be  no  intensity- 

wavelength  window  that  produces  sufficient  symmetrization.  Although 

there  is  experimental  data  with  microvaves^^  that  shows  that  a  broad  bandtdldth 

raises  some  of  the  plasma  instability  thresholds,  there  is  negligible 

data  so  far  with  broad  bandwidth  lasers,  such  as  HF. 

This  work  was  supported  by  the  U.S.  Department  of  Energy  and  the 
Office  of  Naval  Research. 
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